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ABSTRACT 

Rest-frame ultraviolet emission lines offer the exciting possibility to directly image 
the gas around high-redshift galaxies with upcoming optical instruments. We use a 
suite of large, hydrodynamical simulations to predict the nature and detectability of 
emission lines from the intergalactic medium at 2 ^ z ^ 5. The brightest emission 
comes from Hi Lya (1216 A) and the strongest metal line, Cm (977), is about an 
order of magnitude fainter, although Hi Lya may be fainter if the gas is self-shielded 
to the UV background or if dust is important. The highest surface brightness regions 
for CIV (1548,1551), Sim (1207), SilV (1393,1403) and Ovi (1032,1038) are fainter 
than the brightest Cm by factors of a few. The N V (1239,1243) and Ne VIII (770,780) 
lines, as well as Hell Ha (1640), are substantially weaker but their maximum surface 
brightnesses still exceed 10 2 photon s _1 cm~ 2 sr" 1 at z = 2 (for 2" pixels). Lower ion- 
isation lines typically arise in denser and colder gas that produces clumpier emission. 
The brightest H I Lya emission arises exclusively in highly overdense gas, but the high- 
est surface brightness emission from high-ionisation metal lines traces a much wider 
range of overdensities. Bright metal-line emission traces gas with temperatures close 
to the peak of the corresponding emissivity curve. While H I Lya, He II Ha, C III, Si III, 
and Si IV are excellent probes of cold accretion flows and the colder parts of outflows, 
ClV, Nv, OVI, and NeVIII are powerful tracers of the diffuse warm-hot intergalactic 
medium and galactic winds. A comparison of results from simulations with varying 
physical prescriptions demonstrates that the predictions for the brighter metal-line 
emission are robust to within factors of a few. Several rest-frame UV emission lines 
from the high-redshift IGM will become detectable in the near future, possibly starting 
with the Cosmic Web Imager, which is already operating on Palomar. The Multi Unit 
Spectroscopic Explorer, which will be commissioned in 2012 on the Very Large Tele- 
scope, and the proposed Keck Cosmic Web Imager have the potential to revolutionise 
studies of the interactions between high-redshift galaxies and their environment. 

Key words: intergalactic medium - diffuse radiation - radiation mechanisms: ther- 
mal - galaxies: formation - cosmology: theory 



1 INTRODUCTION 

Investigating the spatial distribution of galaxies in the 
sky has taught us much about how the Universe evolves 
and structures form. However, galaxies contain only a few 
percent of the total bar yonic mass in the Universe (e.g. 
iFukugita fc Pee bles 2004j) and their spatial distribution is 
stron gly biased with re s pect to the dark mat ter distribution 
(e.g. lAllen et all 120051 ; blanton et~ai] l200fj ). Mapping the 
distribution of the intergalactic medium (IGM) could pro- 
vide one of the most powerful tools in modern astrophysical 
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cosmology as it contains most of the baryons and is only 
weakly biased with respect to the dark matter. 

While absorption lines in the spectra of background 
quasars provide 1-dimensional information along the line 
of sight, emission lines map the full 3-dimensional distri- 
bution of the gas, thus providing a more realistic picture 
of the physical state of the gas and of its evolution with 
time. Since the IGM is the main reservoir of fresh material 
to fuel star fo rmation in galaxies and is in return affected by 
feedback (e.g.lTheuns et al.ll20o3; iBertone. Stoehr fc Whitel 
20051: IBertone fc White! l200rj; lOppenheimer fc Pavel 12008 
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properties of cosmic gas ultimately implies improving our 
understanding of the process of galaxy formation. 

Predictions for the emission from the low-redshift 
IGM in the rest-frame ultraviolet (UV) have been made 

i 1 r —H I — 

for b oth Hi Lya (Furlanetto ct al. 2003; Furlanctto et al 
2005) and metal lines (jFurlanetto et al. 2004; Be rtone et al 
2010bl . Paper II hereafter). At z > 1.5, the emis- 
sion is redshifted in to the optical band and acces- 
sible to ground-based telescopes, but predictions have 
only been made for Hi Lya and Hen emission 
(e.g. iHoean fc Weymannl Il987l: iGould fc Weinberg! Il996l: 
Fardal et al.l l200ll; IFurlanetto et all 120051; ICantalupo et al. 



20051 ; lYang et al.ll2006l; Diikstra & Loebl 2009; iGoerdt et al 
201ol; iKollmeier et alJ 201Cl: iFaucher-Gie uere et alJl2010h ~ 
ICantalupo et alj (|2005l ) and iKollmeier et alJ |2010l ) in- 
vestigated fluorescent Lya emission induced by the UV 
background and by the ionising radiation of bright quasars. 
IKollmeier et alj l|2010l) claim that while the detection of 
the glow from the diffuse IGM is beyond the sensitivity 
of current telescopes, quasar-induced fluorescence in high- 
density r egions should be detectable around bright quasars 
(see also iFrancis fc Bland-Hawthornll2004l ; ICantalupo et al] 
2005). Indeed, such emission may already have bee n de- 
tected (IFrancis fe McDonnell! 120061 ; I Adelberger et al.l l200 ' 



ICantalupo. Lilly, fc Porcianill2007f ). IFurlanetto et all (| 200 



discussed the relative importance of fluorescence and colli- 
sional excitation and claim, however, that away from bright 
sources the latter is the dominant channel for the production 
of bright Lya emission. lYang et~al ] (|2006h predict that He n 
emission at z as 3 is detectable by current facilities, and 
that the Hen 1640 A line in particular, could be a useful 
tool to investigate the cooling mecha nisms in the so - called 
"Ly a blobs" . Recent observat i ons by IScarlata et al.l |2009) 
and IPrescott. Dev. fc Jannuzil |2009l ) may have confirmed 
this prediction. 

As of today, Lya blobs remain the main form of diffuse 
emission detected a t 
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high redshift ( e.g. ISteidel et al 
iBower et al.l |2004| ; IWeidinger et al 



2000: 



20091 ; IWeiimans et al 



20071; 
120101 ) 



Rauch et al 



2006 



2008 



Saito et al 



Yang et al 



controversial, 
be due to ( 



The Lya 



but their nature is still 
emission in the blobs might 
infalling on to 



radiation from 



haloes (e.g. 


Fardal et all l200ll; 


Furlanetto et al.l l2005t 


Nilsson ct al. 


2006; Diikstra et al. 


20061; Smith & Jarvid 


20071; Diikstra & Loebl 20091: Goerc 


It et alj|201(t), but it is 



difficult to rule out other mechanisms su ch as galactic 
wind s driven by supernova explosions (e .g. lOhvama et al.l 
120031 ; iMori et al.ll2004l ; iGeach et al]|2009l ) and particularly 
recombination radiation powered by internal source s (e.g. 
IFurlanetto et al". I 120051 ; iFaucher-Giguere et all |2010| ). The 
detection (or non-detection) of metal-line emission from the 
same gas phase that produces Lya radiation would help re- 
veal the nature of the blobs. 

In the next few years, a number of optical integral 
field unit spectrographs will come online with the spe- 
cial purpose of detecting diffuse emission fr om the IGM at 
z > 1.5. The Cosmic Web Im ager (CWI, I Rahman et al.l 
120061 : iMatuszewski et"al] l20ld ) started operating on the 
200 inch telescope at Palomar in 2009, while the M ulti Unit 
Spectroscopic Explorer (MUSE. iBacon et al.l l2009h for the 
Very Large Telescope (VLT) and the Keck Cosmic Web Im- 



Table 1. Adopted solar abundanc es, f rom I Al lende P rieto et all 
(200lh, lAl"iende Prieto et al. I J2002F) and lHolweged fcOOll ). 

Element n^/nu Element n,i/nn 



H 

He 

C 

K 

O 

Nc 



1 

0.1 

2.46X10" 4 
8.51xl0~ 5 
4.90xl0~ 4 
l.OOxlO" 4 



Mg 

Si 

S 

Ca 
Fe 



3.47xl0~ 5 
3.47xl0~ 5 
1.86xl(T 5 
2.29xl0~ 6 
2.82xl0~ 5 



ager (KCWI, iMartin et alj 12009] ) will begin operations in 
2012/2013. These will hopefully be follo wed by the propose d 
Antarctic Cosmic Web Imager (ACWI. iMoore et alj [20081 ). 
In addition, balloon experiments in the near UV, such as 
the Fai nt Intergalactic Re dshifted Emission Balloon (FIRE- 
BALL, iTuttle et aLl l2008). are set to explore the Universe 
at lower redshifts. 

These instruments have the potential to revolutionise 
our understanding of galaxy formation. The ability to map 
the surface brightness and velocity profiles of the circum- 
galactic gas in one or more emission lines would provide 
invaluable information about two crucial, but poorly un- 
derstood processes: gas accretion onto galaxies and galactic 
winds driven by feedback from star formation and/or accret- 
ing supermassive black holes. 

In this paper we investigate the rest-frame UV emission 
from the IGM at redshifts 2 < z < 5, and explore the possi- 
bility of detecting it with ground-based, optical telescopes. 
We will focus on metal-line emission, but will provide some 
estimates for Lya as a point of reference. Using a subset 
of simulations fr om the Overwhelm ingly Large Simulations 
project (OWLS. ISchave et alJlioiol ). we will show that up- 
coming instruments such as MUSE and KCWI have the po- 
tential to detect diffuse emission from multiple transitions 
arising in the gas around galaxies. 

This paper is organised as follows. Section [2] briefly de- 
scribes the OWLS runs and our method to calculate and 
map the gas emission. Section|3]presents detailed results for 
emission in the reference (REF, hereafter) model. In Section 
[3] we compare the predictions of the reference model with 
results from a selection of other OWLS runs with differ- 
ent implementations of a number of physical prescriptions. 
Finally, we draw our conclusions in Section [5] and present 
convergence tests in the Appendix. 

Throug hout this pap e r we assume a WMAP3 ACDM 
cosmology l|Spergel et al.l 120071 ) with parameters fi m = 
0.238, fi b = 0.0418, fi A = 0.762, n = 0.951, and ag = 0.74. 
The Hubble constant is parametrised as Ho = 100 h' 1 km 
s _1 Mpc -1 , with h — 0.73. For reference, the cosmic mean 
density corresponds to a hydrogen number density of nu = 
5.1 X 10~ 6 cm~ 3 [(l + z)/3] 3 . We adopt a solar abundance of 
Zq — 0.0127, corresponding to the val ue obtained using th e 
default abundance set of CLOUDY l|Ferland et alj|l998h . 
This abu ndance set, listed in Table [ H combines the abun- 
dance s of I Allende Prieto etldl ^OOlf T lAllende Prieto et all 
(|2002h and lHolwegerl (|200ll ) and assumes riHe/fiH = 0.1. 
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2 NUMERICAL METHODS 

In this Section we discuss the numerical set up of this work. 
In particular, in Se ction | 2. II we briefly describe the OWLS 
simulation project l|Schave et alj|2010h and in Sections 12.2 
and 12.31 the emissivity tables and the method to calculate 
and map the gas emission, respectively. 



2.1 The simulations 

The OWLS project includes more than fifty different cos- 
mological, hydrodynamical simulations with varying mass 
resolution, box sizes, and physical prescriptions for star for- 
mation, galactic winds driven by supernovae (SNe), feedback 
from active galactic nuclei (AGN), chemical evolution and 
gas cooling. All simulations are produced with a significantly 
modified version of t he paralle l PMT ree-SPH code GADGET 
III, last described in ISpringell (|2005T l. In this work we use 
a subset of simulations with a box size of 25 ft -1 comoving 
Mpc that contain 2 x 512 3 particles with initial gas particle 
masses of 1.4 x 10 6 ft -1 Mq and dark matter particle masses 
of 6.3 x 10 6 hT 1 Mq. The evolution of the dark matter and 
gas density distribution is followed from z = 127 down to 
z = 2. Comoving gravitational softenings are set to 1/25 of 
the mean comoving inter-particle separation, but are limited 
to a maximum physical scale of 0.5 kpc, which is reached 
at z = 2.91. The box size and the mass resolution of these 
simulations make them ideal for studying gas on relatively 
large scales, but with sufficiently high resolution to resolve 
galaxies at z ^ 2. In addition to gravitational and hydro- 
dynamical forces, the simulations include prescriptions for 
radiative cooling, star formation, chemodynamics, SN and, 
for one simulation, AGN feedback. Below we will summarise 
the sub-grid physics implemented in our reference simula- 
tion. In Section 2] we discuss a number of other models that 
incorporate variations of the physical modules, implemented 
one at a time. 

Radiative co oling is implemente d follo wing 
IWiersma et al. | (|2009al r1. As iBertone et al. | l|2010al . Pa- 
per I hereafter) and Paper II showed, accurate cooling 
rates are essential to determine the intensity of the gas 
emission. The net radiative cooling rates are computed for 
11 elements (hydrogen, helium, carbon, nitrogen, oxygen, 
neon, magnesium, silicon, sulphur, calcium, and iron) under 
the assumption of an optically thin gas in (photo-)ionisation 
equilibrium, in the presen ce of the cosmic m i crowav e back- 
ground (CMB) and the lHaardt fc Madaul |200ll 'l model 
for the ionising background radiation from quasars and 
galaxies. The contribution of each individual element to the 
cooling rate is calculated by interpolating pre-computed 
tables (created with the publicly avail able photo-ionization 
package CLOUDY, last described by IFerland et all Il998l ) 
as a function of density, temperature, and redshift. 

Sta r formation is implemen t ed fo llowing the prescrip- 
tions of ISchave fc Dalla Vecch ia (2008), according to which 
gas denser than 7iH,crit 0.1 cm -3 (i.e. the interstellar 
medium) obeys an effective equation of state of the form 
P oc p 7ctt , where P is the gas pressure, p its density and 



1 Using their equation (3) rather than (4) and CLOUDY version 
05.07 rather than 07.02. 



7efi = 4/3. Gas on the effective equation of state is allowed to 
form stars at a pressure-dependen t rate that repro duces the 
observed Kennicutt-Schmidt law (|Kennicutt 19981). Chemi- 
cal evo lution is implemented as described in Wiersma et al.l 
(2009bJ). The simulation tracks the evolution of the abun- 
dances of the 11 elements that contribute to the radiative 
cooling and follows their delayed release by massive stars 
(Type II SNe and stellar winds) and intermediate mass stars 
(Type la SNe and asymptotic giant branch stars ), assum- 
ing a Chabrier stellar initial mass function (IMF) |Chabriei1 
2003) spanning the range 0.1 - 100 Mq. 

SN fe edback is injected in kinetic f orm following the 
method of iDalla Vecchia fc Schavel (|2008T ). In this prescrip- 
tion, type II SNe inject kinetic energy locally into the 
gaseous neighbours of newly formed star particles. The feed- 
back prescription includes two parameters that describe, re- 
spectively, the amount of mass injected into winds per unit 
stellar mass formed, 77, and the initial wind velocity, v v . The 
REF model uses r\ = 2 and v„ = 600 km s , which cor- 
responds to about 40 percent of the total available kinetic 
energy from SNe, with the rest assumed to be lost radia- 
tively. The reference model does not include a prescription 
for AGN feedback, but we will consider a model with AGN 
in Section [4] 



2.2 The emissivity tables 

The gas emission is calculated following the method of Pa- 
per I and we refer the reader to that work for the details of 
the procedures, which we summarise briefly in the following. 

Tables of the line emissivity as a function of temper- 
ature, density and redshift were crea ted with CLOUDY 
(version c07.02.02; IFerland et al.l fl998). under the assump- 
tions of solar abundances, optically thin gas, and (photo- 
)ionisation equilibrium i n the presence of the CMB and the 
lHaardt fc Madaul <|200lf ) model for the evolving UV/X-ray 
background radiation from galaxies and quasars. The ta- 
bles include a total of about 2000 emission lines for the 
11 elements explicitly tracked by OWLS. The temperature 
is sampled in bins of ALog 10 T = 0.05 in the range 10 2 K 
< T < 10 8 ' 5 K and the hydrogen number density in bins of 
ALog 10 nH = 0.2 in the range 10 -8 cm -3 < nrr < 10 cm -3 . 

Fig. [T] shows the emissivity of the lines used in this 
paper as a function of temperature at 2 = 2, for nn — 10 -3 
cm~ 3 . The dependence of the line emissivity on density is 



illustrated in Fig. 1 of Paper II. At low temperature and in 
the low-density regime, the contribution of photo-ionisation 
by the UV background on the ionisation balance becomes 
important and may exceed the contribution from collisional 
ionisation. 

The assumption of ionisation equilibrium in the calcu- 
lation of the cooling rates and line emissivities is justified for 
regions that are predominantly photo-ionised and for dense 
gas in the centres of clusters, but might not be entirely 
robust for the warm-hot intergalactic medium (WHIM) 
(lYoshida et al.ll2005l;lYoshikawa fc Sasakill2006l;ICen fc Fand 



20061; I Gnat fc Sternberd 120071 ; IGnat fc Sternberd 120091 ; 



Bertone et al. 20081 forareview). However, since these stud- 



ies ignored photo-ionisation, they may have overestimated 
the importance of non-equilibrium ionisation. 

The assumption that the gas is optically thin to ionising 
radiation may break down for very high densities. The metal 
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Figure 1. The emissivities of a sample of emission lines as a function of temperature, for solar element abundances, hydrogen number 
density njj = 10 — 3 cm -3 and 2 = 2. See Table[3]for the corresponding rest wavelengths. The tail of the curves at the lowest temperatures 
is due to the effect of photo-ionisation on the ionisation balance, which is density-dependent. 



lines are predominantly collisionally ionised at the temper- 
atures corresponding to the brightest emission. This follows 
from the fact that the emissivity curves are very strongly 
peaked and that the peaks occur at temperatures ^> 10 4 K, 
as can be seen from Fig. \T\ Indeed, we will show in Sec- 
tion [3]4] that while the bright metal- line emission occurs in 
gas with temperatures close to the peak of the emissivity 
curve, the density can vary by orders of magnitudes and it 
is not necessarily the densest gas that dominates. For hydro- 
gen and helium, however, the brightest emission does tend 
to occur in the densest gas with T <C 10 J K. We therefore 
expect that self-shielding may be important for hydrogen 
and helium, but probably not for the metal lines. 

As discussed in detail by iFurlanetto et all (120051 ) and 
iFaucher-Giguere et al.l (|2010h . self-shielding from ionising 
radiation may be very important for H i Lya. While a higher 
neutral fraction may boost the Lya emissivity, the absence 
of photo-heating may cause the temperature of the gas to 
drop to values that are too low to excite the electrons col- 
lisionally, thereby strongly suppressing the emissivity. As a 
consequence, it is possible that the Hi Lya emission from 
the densest part of the IGM could be substantially higher 
or lower than we predict. Another radiative transfer effect 
that we will ignore is line scattering, which is important 
for H i Lya and can for example result in enhanced dust 
extinction. On the other hand, scattering of recombination 
radiation from Hn regions, which we also do not consider, 
may boost the Lya flux from the surrounding circumgalactic 
medium. 

Because of the potential i mportance of radi a tive t rans- 
fer for Lya emission, we follow IFurlanetto et al l (2005) and 
consider not only the optically thin case, but also provide 
predictions after setting the Lya emissivity of all gas with 
nu > 10 -3 cm -3 to zero. This second case is highly conser- 
vative, both because we do not expect all Lya photons from 
higher density gas to be destroyed by dust (after all, com- 
pact Lya emitters are known to exist) and because the gas is 



likely to be fully self-shielded only for higher densities, even 
if w e consider only io nisation by the background radiation 
(e.g. lAltav et aUboill ). 

We note that metal-line emission may in principle also 
be affected by self-shielding, although Fig. [T] (see also Paper 
II) suggests that the bright emission is dominated by colli- 
sionally ionised gas and therefore insensitive to the amount 
of self-shielding. Finally, including local sources could en- 
hance t he de gree of ionization and the temperature (e.g. 
Scha vd I2006T). p articularly in the case of X-ray emission 
l |Cantalupdl2010l ). 

2.3 The gas emission 

In this Section we briefly describe the procedures to calculate 
the gas metal-line emission and to create emission maps. 
These follow the prescriptions of Paper I and we refer the 
interested reader to that work for more details. 

The luminosities of gas particles are calculated by inter- 
polating the emissivity tables over hydrogen number density 
and temperature at the desired redshift, and by rescaling 
the emission to the abundance of the corresponding element. 
The procedure is repeated for each emission line. We use the 
"smoothed" element abundances defined in W iersma et al.l 
(2009b), which alleviates the effect of the lack of metal mix- 
ing inherent to SPH. The emission is calculated only for gas 
particles with density nu ^0.1 cm~ 3 , which thus excludes 
all particles on to which we impose an effective equation of 
state. We exclude these particles, which constitute the multi- 
phase interstellar medium in the simulations, both because 
we lack the resolution and physics required to model this 
dense gas and because we are interested in emission from the 
IGM. Our predictions for the surface brightness are there- 
fore conservative as in reality the brightest emission may 
come from denser gas. 

The luminosity in emission line I, in units of erg s _1 , for 
each gas particle i as a function of its gas mass m gas ,i, density 
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Table 2. Map properties. The pixel size is 2". The angular size 
of the field and the velocity width of the slice both correspond to 
a comoving length of 25 h~ 1 Mpc. 



Redshift 


2 


3 


4 


5 


Number of pixels 


671 2 


542 2 


478 2 


439 2 


Pixel size (comov. h^ 1 kpc) 


37.2 


46.1 


52.3 


56.9 


Angular size (arcmin) 


22.4 


18.1 


15.9 


14.6 


Slice thickness (km s — 1 ) 


2234 


2500 


2762 


3010 



pi, hydrogen number density nn,i and element abundance 
X 7t i is computed as 



Li,i (z) = £i, (*,Ti,riH,<) ■ 



x, 



y© 



(1) 



where ei,q (z,T,nn) is the emissivity of the line I, in units 
of erg cm -3 s _1 and for solar abundances, bi-linearly inter- 
polated (in logarithmic space) from the CLOUDY tables as 
a function of the particle temperature Log 10 Ti and hydro- 
gen number density Log 10 nH,i at the desired redshift z. Xy,i 
is the mass fraction of element y and X q its solar value. 
The corresponding particle flux Fij, in units of photon s _1 
cm~ 2 , is given by 



Fit = 



Lij A; 



(l + z) : 



(2) 



with Dl the luminosity distance, hp the Planck constant, 
c the speed of light and A; the rest-frame wavelength of 
emission line I. 

Finally, to create emission maps we project the fluxes 
onto a 2-dimensional grid whose number of pixels is deter- 
mined by the desired angular resolution and by the redshift 
we want to "observe" . The dependence of the results on the 
angular resolution is illustrated in Appendix [X] The pro- 
jection is done using a flux- conserving SPH interpolation 
scheme. The gas surface brightness is computed by dividing 
the flux by the solid angle Q that subtends a pixel, that is 
Sb,i = Fi/Sl. 

Unless stated otherwise, all emission maps are cre- 
ated using a pixel size of 2" x 2" , approximately consis- 
tent with the angular resolution of CWI (2.5" x 1"; 
iMatuszewski et"aj] |201fJl. and a slice thickness of 25 h~ x 
comoving Mpc, equal to the size of the simulation box. In 
Appendix lB3l we show that for sufficiently large fluxes, which 
include the regime of interest here, the probability distribu- 
tion function (PDF) of the pixel flux is proportional to the 
slice thickness (although we do find that projection effects 
slightly enhance the flux for the case of C iv) . For fluxes that 
are potentially detectable the PDF can thus easily be scaled 
to a different slice thickness. For simplicity, we have com- 
puted the PDFs using the entire simulation box. At z = 2, 
a slice thickness of 25 h~ x Mpc corresponds to a spectral 
resolution Av = 2.2 x 10 3 km s~ , which is much worse than 
proposed for upcoming instruments. Even reducing the slice 
thickness to 1 h^ 1 Mpc, as we do in Fig. IB3I the Hubble 
velocity difference across the slice is larger than the reso- 
lution of CWI. The number of pixels in each map depends 
on redshift. An angular resolution of 2" requires a grid with 
671 2 pixels, which corresponds to a comoving size of 37 h~ x 
kpc at z = 2, or 74 times larger than the gravitational force 



resolution used in our simulations. For our cosmology, a re- 
gion with comoving size 25 h~ Mpc corresponds to a field 
of view of about 22' at z = 2. The same quantities are listed 
for other redshifts in Table [2] 



3 RESULTS 

In this Section we present results for the rest-frame UV emis- 
sion from the IGM in the reference model. In particular, 
we investigate how the intensity of the emission lines varies 
with redshift in Section 13.11 and we discuss the detectabil- 
ity in Section 13.21 We study line ratios in Section 13.31 and 
we investigate how the emission depends on the properties 
of the gas in Section 13.41 In Appendix [B] we demonstrate 
that the results are converged with respect to the numerical 
resolution and the size of the simulation box, although it 
is possible that a larger box size would have allowed us to 
sample the PDF to higher fluxes in the regime where it is 
steeply declining. 



3.1 Line emission 

In this Section we investigate the strengths of the emission 
lines listed in Table [3] and their evolution from z = 5 to 2. 
The lines have been selected on the basis of their strength or 
observability in the optical band. The minimum and max- 
imum redshifts for which the various lines fall within the 
wavelength ranges of MUSE and KCWI are listed in Ta- 
bled 

Among the promising metal lines are those from 
lithium-like ions, i.e. Civ (1548,1551), Nv (1238,1243), Ovi 
(1032,1038) and Nevm (770,780), as well as Cm (977), Sim 
(1207), and the sodium-like ion Siiv (1394,1403). We will 
also present results for two important cooling lines of hy- 
drogen and helium: Hi Lya (1216) and Hen Ha (1640). 
Several other rest-frame UV lines, e.g. Hen Lya (304), Niv 
(765), Oiv (549), and Ov (630), are usually about an order 
of magnitude stronger than lines from more ionised atoms. 
However, we do not consider them here because they only 
shift into the optical at very high redshifts. 

The lines that are singlets may be hard to identify un- 
less they are cross-correlated with other lines and/or galax- 
ies. However, these lines are interesting for a number of rea- 
sons. Firstly, they provide a sensitive diagnostic for gas at 
lower temperatures than the standard lines from lithium-like 
atoms. Secondly, the simultaneous detection of lines from 
different ions provides stronger constraints on the ionisa- 
tion state of the emitting gas, and therefore on its density 
and temperature. Finally, we will show in the following that, 
in the absence of strong foreground absorption, the singlets 
C in and Si in are stronger than the corresponding doublets 
C iv and Si iv. 

Figs. [2] and [3] show maps of the surface brightness of all 
lines in a high-density region centered on a forming group. 
The line intensity is mapped at z = 2, 3, 4 and 5, for a pixel 
size of 2". In case of doublets, such as for example Civ and 
Ovi, we have only considered the strongest component. The 
intensity of the weaker line in the doublet can be obtained 
by dividing the intensity of the stronger one by a factor of 
2. Note that not all lines fall within the optical range at 
all redshifts shown. For comparison, maps of the density, 
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Table 3. Summary of the emission lines considered in this work. The first column lists the ion names. The second column lists the 
rest wavelength A of the main atomic transition and the third the wavelength A2 of the weaker line in doublets (or Ly/3 for the case 
of Hi). Column 4 lists the type of transition for the case of Hi and Hell and the atomic state of the ion for the others. The last four 
columns show the redshift range for which the transition falls within the wavelength range covered by MUSE (4650-9300A) and KCWI 
(3500-10,000A). 



MUSE KCWI 



Line 


A (A) 


A 2 (A) 


Atomic state/Transition 


^min 


^max 




z max 


Hi 


1215.7 


1025.7 


Lya 


2.8 


6.6 


1.9 


7.2 


He 11 


1640.0 




Ha 


1.8 


1.7 


1.1 


5.1 


Cm 


977.0 




Be 


3.8 


8.5 


2.6 


9.2 


Civ 


1548.2 


1550.8 


Li 


2.0 


5.0 


1.3 


5.5 


Nv 


1238.2 


1242.8 


Li 


2.8 


6.5 


1.8 


7.1 


Ovi 


1031.9 


1037.6 


Li 


3.5 


8.0 


2.4 


8.7 


Ne viii 


770.4 


780.3 


Li 


5.0 


11.1 


3.5 


12.0 


Sim 


1206.5 




Mg 


2.9 


6.7 


1.9 


7.3 


Si iv 


1393.8 


1402.8 


Na 


2.3 


5.7 


1.5 


6.2 



temperature and metallicity of the gas in the same region 
are shown in Fig. [3] at z = 2 for a pixel size of 0.5". For 
reference, 

1 r\-19 -1 -2 -2 ^ 

10 erg s cm arcsec ~ 

8.57 x 10 2 ( A ° 5 ^ ( i±-^ photon s -1 cnr 2 sr -1 (3) 

VioooAy v 4 / 

where Ao the rest-frame wavelength. 

We showed in Paper II that at low redshift the UV emis- 
sion comes from dense gas in filaments and in the haloes of 
galaxies, but not from the centres of large groups, because 
these contain gas heated to temperatures in excess of 10 6 K 
and emit mostly in the X-ray band. At high redshift, the 
situation is different. At z > 2 massive structures such as 
clusters and large groups have not yet fully assembled. More- 
over, fixed overdensities correspond to higher proper densi- 
ties and thus higher cooling rates. Cold streams are therefore 
more prominent at higher redshift. Most of the UV emission 
comes from gas inside or in the vicinity of collapsing objects, 
such as the haloes of galaxies and the gas that starts to ac- 
cumulate in filaments and groups. Little or no metal-line 
emission is produced by the lowest density IGM. 

Fig. [5] shows the surface brightness probability distri- 
bution function (PDF) of emission lines produced by the 
IGM (i.e. n H < 10" 1 cm -3 ) at z = 2. The left panel shows 
results for Hi Lya, Hi LyaSS (i.e. Hi Lya from gas with 
riH < 10" 3 cm" 3 , as would be appropriate if no Lya pho- 
tons were produced in or escaped from self-shielded gas), 
Hen Ha, Cm, and Civ and the right panel for Nv, Ovi, 
Nevm, Sim and Siiv. In all panels the plot in the main 
window shows the high flux tail of the PDF, while the inset 
shows the full distribution. For numerical convenience, the 
minimum value of the surface brightness has been fixed to 
S B = 10" 15 photon s" 1 cm 2 sr 1 . This is reflected in the 
high values of the PDFs of metal lines for the lowest values 
in the inset plots. Note that these PDFs, as well as the ones 
in other figures, assume 2" pixels and that, different from 
Figs [2] - |4l they have been computed from maps covering 
the entire simulation box. The dependence on the pixel size 
is illustrated in Fig. I All for a subset of the lines. 

Fig. [6] shows the evolution of the flux PDF for individ- 
ual emission lines, comparing the results at z = 2, 3, 4 and 5. 
Figs.[S]and[S]demonstrate that if we assume the gas to be op- 



tically thin, then the strongest rest-frame UV emission line 
among those shown here is Hi Lya, which at z — 2 has an 
intensity that is about an order of magnitude greater than 
the brightest C in, which is itself a factor of a few brighter 
than the strongest C iv. The difference becomes even larger 
at higher redshifts, where the emission from metals is sup- 
pressed because of the lower gas metallicity. If, instead, we 
conservatively assume that no Hi Lya comes from gas with 
nu > 10" 3 cm" 3 (case Hi LyaSS), then the bright-end of 
the flux PDF is cut off at ~ 10 2 photon s" 1 cm -2 sr -1 and 
the strongest emission is dominated by metal lines. From the 
metal lines shown here, the brightest emission is produced 
by C iv, Si iv and O vi among the higher ionisation states 
and by C in and Si m among the lower ionisation states. 

From Fig. [5] we see that C in and C iv trace each other 
relatively well, with C m somewhat stronger and more bi- 
ased to higher density gas. The same is true for Sim and 
Si iv (Fig. [3J), which, however, fall off more rapidly with 
declining density than the carbon lines. The other metal 
lines, and O vi in particular, are distributed more smoothly. 
The maps illustrate clearly that lines from higher ionisa- 
tion states are better tracers of lower density regions like 
filaments. In practice, the emission of all lines outside col- 
lapsed regions is extremely low and well below the direct 
detection limits of planned detectors. 

The intensities of Nv and Nevm lines are lower than 
those of oxygen and carbon by at least an order of magni- 
tude. Nv might be barely strong enough to be detectable 
at z « 2. While Nevm lines are unlikely to be detectable 
by CWI, they might be within reach of KCWI. However, at 
z < 3.5 a near UV instrument is required. 

Fig. [6] shows that the maximum value of the predicted 
flux increases with decreasing redshift for all lines. For exam- 
ple, at z = 5 the maximum O vi intensity is about 10 photon 
s _1 cm -2 sr -1 , while it is about 300 photon s _1 cm -2 sr -1 
at z = 2, i.e. more than an order of magnitude larger. At 
the highest redshift a significant number of pixels in the 
map do not contain any metal-line emission, as the spikes 
in the PDFs at the minimum value 10 -15 photon s _1 cm -2 
sr" 1 demonstrate. However, as time proceeds, the fraction 
of pixels with very low or no metal-line emission strongly 
decreases, indicating the progression of the IGM metal en- 
richment. The weak evolution of the intensity of the Lya ra- 
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pixel = 2 arcsec 




-2-10 1 2 3 

Log 10 S B (photon cm" 2 s"' sr" 1 ) 



Figure 2. The rest-frame UV emission of the IGM (i.e. gas with njj < 10 -1 cm -3 ) at redshift z = 2 (upper panels), 3 (upper-middle 
panels), 4 (lower-middle panels) and 5 (bottom panels). From left to right we show emissions from Hi Lya, Hi LyoSS (i.e. Hi Lyo from 
gas with njj < 10 — 3 cm -3 ), Hell H«, Cm, and Civ. The maps show a zoom of a region of 4 ft -1 comoving Mpc on a side that contains 
a group forming at the intersection of a web of filaments. All maps use a pixel size of 2" , which corresponds to comoving sizes of 37, 46, 
52 and 57 h~ 1 kpc at z = 2, 3, 4 and 5 respectively, and comoving slice thickness 25 h -1 Mpc. The number of pixels shown in each panel 
is 108 2 , 87 2 , 76 2 and 69 2 at z = 2, 3, 4 and 5, respectively. A given colour corresponds to the same surface brightness level in all panels. 
The minimum and maximum values used for the colour scale do not correspond to the actual minimum and maximum. 



diation, with variations w ell within an order of m agnitude, 
is in good agreement with iFurlanetto et al.l |2005l ). 

The shape of the flux PDF varies for different lines and 
changes with time. The PDFs of the O vi and Ne vni lines 
are narrow and peak at > 1CP 5 photon s" 1 cm~ 2 sr -1 at 
all redshifts. For most other lines, and particularly for C III, 
Si in, and Siiv, the peak of the distribution moves to low 
fluxes (1CT 10 - 10~ 5 photon s" 1 cm -2 sr -1 ) and the distri- 
bution itself tends to flatten out as redshift decreases. The 
behaviour of the PDFs of the Civ and Nv lines is some- 



what intermediate between these two cases. This can be ex- 
plained by the fact that the emission lines trace different gas 
phases with different spatial distributions, which ultimately 
are responsible for the shape of the PDF. As we will discuss 
in more detail in the next Sections, Ovi and Nevin lines 
trace warm-hot gas with T ~ 10 6 K, while most other lines 
trace warm gas with T < 10 5 K. The fraction of the vol- 
ume filled by shock-heated gas at T ~ 10 6 K increases with 
time, while the metal-enriched, warm gas with T < 10 5 K 
that produces Sim, Siiv, Cm and Civ emission becomes 
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pixel = 2 arcsec 




■4 ' ■ ■ . . 



-2-10 1 2 

Log 10 S B (photon cm" 2 s"' sr" 1 ) 

Figure 3. Same as Fig. [2] but for Nv, Ovi, Nevm, Si III, and Si IV, from left to right. 



progressively more confined to the haloes of galaxies and 
small groups as the Universe expands. As a consequence, the 
PDFs of these lines flatten and move toward lower fluxes. 
The behaviour of Nv lines, produced by gas at T ~ 10 5 K, 
is intermediate between these two cases. 



3.2 Detectability 

A number of optical instruments will come online within 
the next few years with the special purpose of detect- 
ing d iffuse emission from the IGM at z > 1.5. In 2009 
CWI (|Rahman et al.ll2006l ; |Matuszewski et al.ll2010h started 
operating on the 200 inch telescope at Palomar. MUSE 
l|Bacon et alj|2009h will be com missioned at the VLT in 2012 
and KCWI (|Martin et 11120091 ) will hopefully be installed at 



Keck soon after that. In the more distant future these instru- 
ment s may be followed by the proposed ACWI (|Moore et al.l 
2008), which would operate from Antarctica. In addition, 
the exploration of the low-redshift Universe ha s begun with 
the n ear UV balloon experiment FIREBALL (jTuttle et al.l 
|200S| ). 



The main characteristics of CWI, MUSE, KCWI and 
FIREBALL are listed in Table H FIREBALL only covers a 
narrow wavelength range in the UV, 1970 - 2130A, but has 
a large FoV of 160" x 160" (angular resolution 10"). CWI 
has a 60" x40" FoV sampled by 24 2.5" x40" reflective slits. 
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REF_L025N512 pixel = 0.5 arcsec z = 2 
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Figure 4. Maps of the density (left panel), temperature (middle panel) and metallicity (right panel) of the gas in the REF simulation 
at z = 2. The displayed region is the same as that shown in Figs. [2] and \3\ It shows a zoom of a region of 4 h^ 1 comoving Mpc on the 
side that contains a group forming at the intersection of a web of filaments. All maps use a pixel size of 0.5", which corresponds to a 
comoving size of 9 /i -1 kpc at z = 2, and slice thickness 25 h~ 1 comoving Mpc. 



Log, S B (photon s" 1 cm"* arcsec" 2 ) Lc>9io S B (photon s"' cm" 1 arcsec"*) 

-10 -8 -6 -4 -10 -8 -6 -4 




02460246 

Log 10 S B (photon s"' cm" 2 sr"') Log 10 S B (photon s"' cm" 2 sr"') 

Figure 5. The z = 2 surface brightness PDFs for emission from the 2 = 2 IGM (i.e. gas with njj < 10 — 1 cm -3 ) for our sample of 
rest-frame UV emission lines, also shown in Figs. [2] and [3] and listed in Table [3] The left panel shows results for Hi Lyct, Hi LyaSS, 
Hell Ha, Cm, and Civ, and the right panel for Nv, O VI, Nevm, Sim, and Si IV. The main plotting area shows only the high flux tail 
of the distributions, while the full PDFs are shown in the inset. The pixel size is 2" (which corresponds to a comoving size of 37 h" 1 kpc) 
and the PDFs are computed from maps covering the entire simulation box. In Appendix IB3I wc show that for sufficiently large fluxes 
(> 0.1 photon s — 1 cm -2 sr -1 for Ovi) the PDF is proportional to the slice thickness, which we have taken to be 25 h~ 1 comoving Mpc 
(i.e. the size of the simulation box). 
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Figure 6. As Fig. [5] but for Hi Lya, Hi LyaSS, Hen Ha, Cm, and Civ (upper panels, from left to right), Nv, Ovi, Nevin, Sim, 
and Si IV (bottom panels). Each panel shows the results for z = 2,3,4 and 5. The surface brightness on the ir-axis is in units of photon 

—1 —2 —1 

s cm sr . 



Table 4. The characterist ics of the FIR EBALL balloon ex periment llTuttle et alj|2008h as well as the CWI {Matuszewski et al.|[201(t) , 
KCWI {Martin et al.ll2009h . and MUSE {Bacon et al]|2009l) instruments. 



Instrument 


Mirror diam. 


Wavelenj 


;th range 


Bandwidth 


FoV 


Angular resolution 


Resolving power 




(m) 


(A) 


(A) 


(arcsec 2 ) 


(arcsec) 




FIREBALL 


1 


1970 - 


- 2130 


160 


160 x 160 


10.0 


5000 


CWI 


5 


3800 - 


- 9500 


150 


60 x 40 


2.5 X 1 


> 5000 


KCWI 


10 


3500 - 


10000 


300 - 2500 


20 x (8.4 - 34) 


0.35 - 1.4 


900 - 20000 


MUSE 


8 


4650 - 


- 9300 


4650 


60 x 60 


0.2 


1750 - 3750 



Although the wavelength range is Iargfl the instantaneous 
bandwidth is limited to 150A. The KCWI will have various 
modes. The highest spectral resolution mode, R ~ 20, 000, 
will only be available with a small FoV of 20" x8.4" (angu- 
lar resolution 0.35") and for a narrow, instantaneous band 
pass of 300A. On the other hand, at the coarsest spectral 
resolution, R « 900, the FoV is 20" x34" (angular resolution 
of 1.4") and the band width 2500A. There will be a number 
of intermediate settings as well. While MUSE has a much 
bigger FoV (60" X 60"), a higher angular resolution (0.2", 
note that MUSE will make use of adaptive optics) and a 
larger instantaneous band width (4650 — 9300A), its maxi- 
mum resolving power is limited to R = 3750 and it will not 
cover the blue. 

For a spectrally unresolved emission lineB an 80 hour 
observation with MUSE will reach a limiting surface bright- 

2 Currently CWI's band width is still limited to 4400 - 5600A 
{Matuszewski et aljbuld l 

3 The spectral resolution of MUSE varies with wavelength but 
is typically ~ 10 2 km s — 1 . Using a slice thickness that matches 
the spectral resolution does not change our conclusions for metal 
lines. For H I Lya the situation may, however, be different because 
resonant scattering, which our simulations do not model, may 
strongly broaden the lines. 



ness of 2 X 10 19 erg s 1 cm 2 arcsec 2 for a S/N of five 
and a 1 arcsec feature at A = 5000 A. We expect KCWI 
to reach similar depths. Using equation Q, we see that 
MUSE and KCWI should be able to detect intensities of 
~ 2 x 10 3 photon s" 1 cm" 2 sr -1 from the ground for 1" 
features and a S/N of five. For more extended emission 
fainter levels could be reached, of course all provided that the 
background can be subtracted to the required accuracy. As 
some of the emission will indeed be extended (see Fig. IA1[I . 
we thus expect deep observations with upcoming instru- 
ments to be able to directly detect diffuse emission down 
to ~ 10 2 - 10 3 photon s" 1 cm 2 sr 1 . Stacking galaxies 
may allow fainter emissi on to be detected statistically. In- 
deed, UteidereraLl |201ll ) have already detected diffuse Lya 
emission down to 10 _19 erg s _1 cm -2 arcsec -2 by stacking 
z ~ 2.65 star-forming galaxies. 

Table [3] lists the minimum and maximum redshifts for 
which the various lines that we investigate fall within the 
range of MUSE and KCWI, respectively. The redshift range 
for CWI is only slightly smaller than that of KCWI (see 
Table 13). 

According to our predictions, a number of emission 
lines should be detectable by upcoming instruments. As 
a benchmark, we use the instrum ent sensitivity of 200 
photon s _1 cm -2 sr -1 , reported by iRahman et all (|2006l ) 
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for CWI for 10 a r csec features. We note, however, that 
iMatuszewski et"aD (2010) quote a much worse limiting sen- 
sitivity of 10~ 18 erg s" 1 cm" 2 arcsec -2 for features of 
600 arcsec 2 and a 16 hour integration, which corresponds 
to ~ 10 4 photon s" 1 cm -2 sr" 1 . They note, though, that 
the sensitivity can be improved by over an order of mag- 
nitude if the background subtraction is improved using the 
nod-and-shuffie technique. In any case, 200 photon s" 1 cm -2 
sr" 1 will for example be reachable with MUSE for ~ 10", 
spectroscopically unresolved features (at 5cr in an 80 hour 
integration) . 

Hi Lya emission from regions with p > 10 3 pb should 
be detectable up to at least z — 5 (taking our optically thin 
calculations at face value), with the glow of less dense regions 
(p ~ 10 3 pb) visible only at lower redshifts. The Hen Ha line 
may be visible in dense (but still «h < 10" 1 cm -3 ) gas at the 
lowest redshift we consider (z ~ 2). Civ, Ovi and Siiv will 
all be detectable up to z = 3, while Cm and Sim should 
be detectable up to z ~ 4 or higher. The Nv and Nevm 
lines should also be within the capabilities of MUSE and 
KCWI and with some luck they might even be able to detect 
weak O iv and O v emission at z > 4.5. In practice, however, 
detection of Ovi, Ov and Nevm, which have rest-frame 
wavelengths shorter than the hydrogen Lyman limit (912 A), 
may not be feasible due to the high optical depth provided 
by intervening absorption systems. The same may be true for 
Cm and Sim, which have rest-frame wavelengths of 977.0 
and 1206.5 Aand are thus susceptible to intervening Lya 
absorption. As we noted earlier, statistical detections, e.g. 
by stacking H i Lya emitters, should enable the detection of 
fainter lines. 

Recall that since our focus is the IGM, we excluded all 
gas with densities nu > 0.1 cm -3 . As the emissivity scales 
as density squared, this high-density gas may be more easily 
detected, particularly for instruments with high angular res- 
olution such as MUSE and KCWI. On the other hand, for 
observations with instruments with relatively poor angular 
resolution, such as CWI, it may be difficult to tell whether 
apparently diffuse emission is in fact due to a few unresolved, 
dense knots, thus potentially mistaking emission from the 
ISM of small galaxies for emission from diffuse, intergalactic 
gas. 

It will be necessary to reach fluxes well below the optical 
background light (OBL, hereafter) to detect emission from 
the diffuse IGM. Most of the OBL luminosity is due to zo- 
diacal light, with only relatively minor contributions from 
the extragalactic background light and the diffuse galac- 
tic emission. The intensity of the OBL has been measured 
to be 33.5 x 10" 9 , 105.7 x 10" 9 , and 72.4 x 10" 9 erg s" 1 
cm" 2 sr" 1 A" 1 at 3000, 5500, and 8000 A, respectively 
(I Bernstein! l2007i ; see also iBernstein. Freedman Sz Madord 
2002 and references therein). Assuming a resolving power 
of R = A/AA = 10,000, these correspond to ~ 2 x 10 3 , 
2 x 10 4 and 2 x 10 4 photon s" 1 cm" 2 sr" 1 at 3500, 5500 and 
8000 A, respectively. 

As di scussed in Paper II, the FIREBALL balloon ex- 
periment |Tuttle et alj|2008h provides a good chance to de- 
tect IGM emission in the UV band at lower redshifts and 
shorter wavelengths than CWI. The instrument targets the 
wavelength window 1970-2130 A, where atmospheric absorp- 
tion is minimal. The surface brightness limit is about 2000 
photon s" 1 cm" 2 sr" 1 for an angular resolution of 10". The 



instrument aims to detect H i Lya, C iv and O vi emission at 
2 < 1. According to our predictions, in this energy range it 
might be possible to detect the O v line at 630 A at z ~ 2.2, 
if intervening absorption is not very strong. For this emis- 
sion line at this redshift we predict a maximum intensity 
> 10 3 photon s" 1 cm" 2 sr" 1 , which might be within reach of 
the instrument capabilities. FIREBALL might also be able 
to detect Cm emission at z « 1, whose surface brightness 
should be about an order of magnitude larger than its detec- 
tion limit. Unfortunately, both the O v and the C in lines are 
singlets and might thus be difficult to identify in spectra with 
limited redshift and wavelength coverage. As such, it is likely 
that both lines will simply be treated as unknown contami- 
nants. However, cross-correlation with stronger lines, such as 
Hi Lya, may allow a statistical detection of singlets. A pos- 
itive detection would provide proof that C in and O v, and 
lines from beryllium-like atoms in general, are indeed an ef- 
ficient cooling channel and would encourage further studies 
of emission in the far UV band. 

Finally, we caution the reader that the discussion pre- 
sented in this section should only be taken as a rough guide 
to what may be possible with upcoming facilities. More ro- 
bust conclusions would require a proper study of detectabil- 
ity, which would involve analysing virtual observations tar- 
geted to specific instruments, instrumental modes, and ob- 
serving strategies and using realistic analysis techniques. 
While such a study is clearly beyond the scope of the present 
work, the results presented here suggest that such investi- 
gations are worthwhile for several transitions and redshifts. 

3.3 Line ratios 

In this Section we compare the intensity of pairs of lines. 
Results for a subsample of emission lines at z — 2 are pre- 
sented in Fig. [7] for the REF model. The chosen lines are 
representative of all interesting transitions. Fig. [7| focuses 
only on the high surface brightness end of the distribution, 
with Sb ^ 10" 1 photon s" 1 cm" 2 sr" 1 . 

As we will discuss further in the next Sections, the in- 
tensity of the emission lines depends on the gas tempera- 
ture, density and metallicity. At a fixed temperature, suffi- 
ciently high for collisional ionisation to dominate over photo- 
ionisation, the relative intensity of one line with respect to 
another is independent of the density and proportional to 
the ratio of the elemental abundances. Since the shape and 
values of the emissivity curves as a function of temperature 
are characteristic for each line, the line ratios vary line by 
line. This is seen in Fig. [7] The lines of different ions of the 
same element are strongly correlated, but the line ratios of 
ions from different elements show a larger scatter. This is 
because the gas that produces strong emission in one line 
is often different from gas that produces emission in a dif- 
ferent line. Indeed, the scatter is largest for ratios between 
lines whose emissivities peak at very different temperatures. 
Hi Lya emission is almost always stronger than any other 
line (but note that this changes if we ignore Hi Lya from 
self-shielded gas; see e.g. Fig. [5]) , with the exception of a 
few pixels in which carbon and oxygen lines dominate. These 
likely trace highly metal enriched gas in the haloes of galax- 
ies and groups. 

The ratio between the intensity of pairs of emission lines 
ultimately provides information about the ionisation state 
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Figure 7. Comparison of line intensities for a subsample of emission lines, namely Hi Lyo, Hell He*, Cm, Civ, Nv, Ovi, Si III, and 
Si IV. The contours are spaced logarithmically by 1 dex. The distribution of line ratios is shown only for the high surface brightness tail 
of the PDF, i.e. intensities higher than 0.1 photon s _1 cm -2 sr — 1 . The correlation is roughly linear for lines from the same atom (Cm 
and Civ, Sim and Siiv). The scatter in the distributions varies from ratio to ratio, and is largest for lines with emissivities which peak 
at very different temperatures. 



and the temperature of the gas. The intensity of a line is 
highest when the gas has a temperature close to the peak 
temperature of its emissivity curve. A higher value of the 
line ratio indicates that the gas temperature is closer to the 
peak temperature of the line in the numerator. 



3.4 What gas produces the emission? 

In this Section we investigate the physical properties of the 
gas that produces most of the emission. In particular, we 
want to understand what the typical densities and temper- 
atures are of the gas that produces the bulk of the emission 
for each line. This will for example tell us whether the gas is 
bound to galaxies, resides in groups, or traces lower density 
regions like filaments and voids. 
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Figure 9. The emission-weighted gas distributions in the temperature-density plane at z = 2 for the REF simulation. From left to 
right, the panels show the distributions weighted by Hi Lya, Hi LyaSS, Hen Ha, Cm and Civ emission (upper panels), and by Nv, 
Ovi, Nevin, Si III and Si IV emission (lower panels). The horizontal and vertical lines are as in Fig. [H] and the contours are spaced 
logarithmically by 1.5 dex. Carbon and silicon lines trace mostly gas with T < 10 s K, while N V and O VI better trace the WHIM with 
10 5 K< T < 10 6 K. Nevin traces hotter gas, on average, than any other emission line in our sample, with 10 5 K < T < 10 7 K. The 
hydrogen and helium lines trace the gas mass better than the metal lines, although they are biased to cool, dense gas in close proximity 
to galaxies. 
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Figure 8. The distributions of gas (left panel) and metal mass 
(right panel) in the temperature-density plane at z = 2. Results 
are for the REF simulation. The contours are spaced logarith- 
mically by 1.5 dex. The vertical line indicates the threshold at 
nn = 0.1 cm -3 above whi ch we impose an effective eq uation 
of state for star-forming gas {Schavc fc Dalla Vecchiall200sh . The 
horizontal lines define the range of temperatures of WHIM gas, 
10 5 K < T < 10 7 K. 



Fig. [5] shows the mass-weighted (left panel) and the 
metal mass-weighted (right panel) gas distributions in the 
temperature-density plane at z = 2. Most of the mass re- 
sides in intergalactic gas with low density and temperature. 
Compared with the gas mass distribution, the distribution 
of metals in the IGM is biased towards higher densities and 
temperatures. As a consequence, the bulk of the metals do 
not trace the bulk of the mass. Most of the metals are con- 
centrated in two separate areas of the distribution: i) in the 
dense (nH > 0.1 cm -3 ) star-forming gas, i.e. the ISM, for 



which we impose an effective equation of state and which 
we exclude from the analysis, and ii) in warm-hot gas in 
overdense regions with 10 -6 cm -3 < nn < 10 -3 cm -3 . This 
result is qualitati vely similar to the find ings of Paper I for 
z = 0.25 (see also lWiersma et alj[2009bh . 

Fig. [9] shows the emission-weighted gas distributions 
in the temperature-density plane at z = 2. From left to 
right, the panels show the distributions weighted by H i Lya, 
H i Ly«SS, He ii Ha, C in and C iv (upper panels), N v, O vi, 
Nevin, Si in and Siiv (bottom panels). 

A comparison of Figs. [8] and [9] demonstrates that the 
bulk of the emission from metal lines traces neither the IGM, 
nor the metal mass. This is consistent with what we found 
in Papers I and II at lower redshifts for X-ray and UV emis- 
sion lines, respectively. Instead, each metal line traces a spe- 
cial region of the temperature-density plane centered on the 
temperature for which the line emissivity peaks (see Fig.[l|. 
The width of the distribution in temperature corresponds 
roughly to the width of the emissivity curve as a function 
of temperature. In terms of density and metallicity, all the 
emission-weighted distributions peak where the density is 
high, but also where both the metallicity and the density 
are high at the same time. The situation is different for the 
hydrogen and helium lines, which on average trace the gas 
mass well, once the quadratic dependence of the emissivity 
on the density is taken into account. 

There are some notable differences between the different 
emission lines. Carbon and silicon lines trace gas with lower 
temperatures (T < 10 5 K) than other lines. O vi traces dif- 
fuse gas with temperatures in the range 10 5 K < T < 10 6 K. 
N v lines arise from gas with temperatures that fall between 
those of carbon and oxygen lines, while the Nevin lines 
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Figure 10. As Fig. [5] but including only gas in fixed density intervals. From left to right the panels show the PDFs of the pixel surface 
brightness for Hi Lya, Civ, and Ovi. While brightest Hi Lya emission arises exclusively in highly overdense gas. a much wider range 
of overdensities contributes to the brightest C IV and O VI pixels. 



Log t Sr (photon s"' cm" 2 orcsec" 2 ) 




12 3 4 5 6 
Log 10 S B (photon s"' cm" 2 sr"') 



Log 10 S B (photon s"' cm" 2 orcsec" 2 ) 
-10 -9 





-15-10 -5 5 
Log 10 S B 

C 4 

Log,„ T<4 

4<Log 10 T<5 

5<Log 10 T<6 
6<Log 10 T<7 
Log, T>7 



1 2 3 4 5 6 
Log 10 S B (photon s"' cm" 2 sr"') 



Log 10 S B (photon s"' cm" 2 orcsec 




12 3 4 5 6 7 
Log )0 S B (photon s"' cm" 2 sr"') 



Figure 11. As Fig. 1101 but including only gas in fixed temperature intervals. For these metal lines the brightest emission comes from 
gas with temperatures close to the peak of the emissivity curves (see Fig.[TJ. 



trace hotter gas than any other line we consider here, with 
T~ 10 6 K. 

Figs. UOland 1111 show explicitly how the surface bright- 
ness PDF varies with the gas density (Fig. 1 1 Op and tem- 
perature (Fig. 1 11 f) for Hi Lya (left columns), Civ (middle 
columns), and O vi (right columns). We created these figures 
from emission maps that were computed by only including 
particles with the physical properties indicated in the leg- 
ends. While Hi Lya is highly biased towards high-density 
gas, with the highest surface brightness due to gas with 
overdensities p/ (p) > 10 3 , the brightest metal-line emission 
arises from gas with a wide range of overdensities, though 
still > 10. The gas responsible for the strongest C iv and 
Ovi emission has temperatures that are close to the peaks 
of the emissivity curves shown in Fig. [TJ T ~ 10 5 K and 
~ 10 5 ' 3 K, respectively. 



4 VARYING THE PHYSICS 

In this Section we discuss the influence of different physical 
processes on the results. We focus our analysis at z = 2, but 



all our conclusions are valid for results at higher redshifts. 
We consider five different sets of OWLS runs, in which the 
physical modules have been changed one a t a time. The mod- 
els are as follows (see ISchave et al.l 120101 for more detailed 
descriptions) : 



(i) REF: Our reference model, as described in Section [2~T1 
All results presented in the preceding sections were obtained 
from this si mulation. This model uses the st ar formation pre- 
scription of ISchave fc Dalla Vecchial (l200St), the supe rnova- 
driven wind model of iDalla Vecchia fc Schavd (|2008r ) with 
initial wind mass loading r\ — 2 and initial wind velocity 
d w = 600 km s -1 , th e metal-dependent cooling rates of 
W lcrsma et al l J2009al). the chemical evolution model de- 
scribed in lWiersma etaLI l)2009bl 1 and the WMAP3 cosmol- 
ogy. AGN feedback is not included. 

(ii) NOSN.NOZCOOL: As REF, but without supernova 
feedback and with radiative cooling rates that assume pri- 
mordial abundances. The metals produced by stars are thus 
only transported by gas mixing and no energy is transferred 
to the IGM and ISM when SNe explode. 
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Figure 12. As Fig. [2] but for a selection of OWLS models with different physics implementations and for a selection of emission lines at 
z = 2. From top to bottom, the lines shown are H I Lya, C IV and O VI. From left to right, the results are for the REF, NOSN_NOZCOOL, 
NOZCOOL, WVCIRC and A GN models respectively. 



(iii) NOZCOOL: As REF, but with radiative cooling rates 
calculated assuming primordial abundances. Neglecting the 
contributions of metal lines implies slower gas cooling than 
in the REF model. 

(iv) WVCIRC: As REF, but with "momentum- 
driven" galactic win d s fo llowing the prescription of 
lOppenheimer fc Pavel ([2008) , with the difference that, 
as in the REF model, winds are "local" to the star 
form ation event and are fu l ly hy drodynamically coupled 
(see iDalla Vecchia fe Schavel 2008 for a discussion of the 
importance of these effects). The wind velocity and mass 
loading factor are functions of the galaxy velocity dispersion 
a: v„ — 5a and 77 = « w o/cr, with v w o = 150 km s _ . The 
velocity dispersion a = \/2\>circ is estimated using an 
on-the-fly friends-of-friends halo finder. This model was 
motivated by the idea that galactic winds may be driven 
by radiation pressure on dust grains rather than SNe. The 
energy injected into the wind becomes much higher than 
for REF for halo masses exceeding 10 11 — 10 12 Mq. 

(v) AGN: As REF, but wit h the addition of AGN feed- 
back, using the prescription of lBooth fc Schavel (|2009l ) . 



We chose this set of simulations because these are the 
models that produced the largest variations in Papers I and 
II. In these papers, which focused on 2 = 0.25, we found 
that for most models the variations in the results are small 
relative to the REF model, especially for UV lines. In terms 
of the bright-end tail of the surface brightness PDF, the 
largest discrepancies with the REF results were found for 
model NOZCOOL. The inclusion of AGN feedback could 
decrease the X-ray emission at the centres of groups and 
clusters by up to a factor of ~ 100, but barely affected the 
emission from diffuse IGM gas. Paper II also showed that 
varying the feedback parameters that control the intensity 
of the winds (e.g. the wind mass loss rate and velocity) in- 
troduced variations in the UV results of a factor of a few or 
less. 

Fig. [T2]shows maps of the surface brightness of H 1 Lya, 
C iv and O vi lines in the different simulations. The PDFs 
of the surface brightness are shown in Fig. [13] All maps as- 
sume an angular resolution of 2" and show results at z — 2. 
Most models predict comparable maximum surface bright- 
ness values, as can be seen by comparing the high surface 
brightness tails of the PDFs (Sb > 100 photon s _1 cm -2 
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Figure 13. As Fig. \E\ but for a selection of OWLS models that include different physical processes. From left to right, the different 
panels are for z = 2 emission from H I Lya, C IV and O VI, respectively. There is almost no variation in the H I Lya emission predicted 
by the different simulations, although removing emission from gas with njj > 10" 3 cm - 3 , as would be appropriate if no Lya photons 
would be produced in or escape from such regions, would cut off the bight end of the distribution. The differences are largest for low 
surface brightness metal-line emission, but the models tend to converge at the bright end. Metal-line cooling is particularly important 
for Ovi. AGN feedback slightly decreases the peak surface brightness of both Civ and O VI. 



sr -1 ), but the variations can be very large for fainter metal- 
line emission. For H i Lya, on the other hand, the differences 
are always very small. 

In the NOSN.NOZCOOL model the emission is concen- 
trated close to galaxy haloes, while the surface brightness of 
low-density regions is extremely low for metal lines. This 
produces the flattening observed in the PDFs of metal lines 
for Sb < — 100 photon s" 1 cm" 2 sr" 1 and is due to the 
inefficient transport of metals from galaxies to the diffuse 
IGM in the absence of SN feedback. 

Model NOZCOOL predicts somewhat stronger emission 
from metals than the REF model, but comparable H i Lya 
emission. In dense regions, the O vi emission exceeds the pre- 
dictions of the REF model by up to an order of magnitude. 
This indicates that metal-line cooling is important in the 
gas for which the Ovi emission peaks, i.e. highly enriched 
gas with T ~ 10 5 ' 3 K. Note that ignoring metal-line cooling 
is not self-consistent, as it allows the gas to linger too long 
at the temperatures for which the emissivity is high. 

For model WVCIRC the peaks of the metal-line PDFs 
are shifted to higher surface brightness values. In other 
words, compared with model REF, a larger fraction of the 
sky, and a larger fraction of the mass in the Universe, shines 
bright in the UV. The outflows in model WVCIRC axe more 
efficient at transporting metals into the IGM, mainly as a 
result of the higher mass loading in low-mass galaxies. As 
such, the different spatial distribution of the UV emission 
we see in Fig. [12] between the REF and WVCIRC models is 
mostly due to differences in the metallicity of the gas, and 
to a lesser extent to differences in the temperature and den- 
sity. In principle, the detection of diffuse metal-line emission 
could thus help decide which one of the descriptions of the 
physics of galactic winds is more "realistic". A larger spa- 
tial extent and a more filamentary structure of the emission 



on large scales would weigh in favour of highly mass-loaded 
galactic winds. 

Finally, the predictions of model AGN are compara- 
ble to those of the REF model in regions with low surface 
brightness and low density. However, at the bright tail of the 
distribution model AGN predicts surface brightnesses that 
are up to 0.5 dex lower. This is a consequence of the fact 
that AGN feedback decreases the density and metallicity of 
the gas near massive galaxies. 

The relatively small differences between the metal-line 
emission predicted by models that include very different 
physical processes gives us confidence that our predictions 
for the brighter emission are robust to factors of a few. On 
the other hand, the fact that the differences in the predic- 
tions for the diffuse metal-line emission are not negligible, 
implies that observations of this faint glow have the poten- 
tial to discriminate between different models for feedback 
from star formation and AGN. 

Finally, we recall that while the different simulations 
predict very similar Hi Lya emission, the predictions for 
this line are uncertain due to the importance of self-shielding 
from ionising radiation, resonant scattering, and absorp- 
tion by dust, which are all ignored here. The potential im- 
pact is illustrated by the blue dashed curve labelled 'REF 
SS' in the left panel of Fig. 1131 which indicates the sur- 
face brightness PDF if we ignore all emission from gas with 
nu > 10~ 3 cm~ 3 , an assumption that is very conservative 
for our purposes. 



5 CONCLUSIONS 

We have investigated the nature and detectability of red- 
shifted rest-frame UV line emission from the IGM (i.e. 
nu < 10" 1 cm" 3 ) at 2 < z < 5. While more strongly bi- 
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ased towards high-density gas than absorption lines, emis- 
sion lines offer the possibility to directly map the distribu- 
tion of the gas a nd metals in three dim ensions. Following 
the procedures of iBertone et all l|2010al ). we used a set of 
large, cosmological, hydrodynamical sim ulations taken from 
the OWLS project (ISchave et al.l|20 10) to create simulated 
maps of a set of strong rest-frame UV emission lines. We 
focused on Cm, Civ, Nv, Ovi, Nevin, Sim and Siiv lines 
(see Table [3}, which are the strongest metal lines that trace 
the diffuse IGM and that can potentially be detected by 
current and upcoming optical instruments. For reference, 
we also showed results for H i Lya and He n Ha and for the 
former we illustrated the potential impact of self-shielding 
and dust extinction by considering only the flux emitted by 
gas with riH < lCP 3 cm -3 , which provides a conservative 
estimate of the minimum signal that can be expected. 
Our main results can be summarised as follows: 

(i) Depending on the line and the redshift, we predict 
a maximum surface brightnesses of order 10 1 — 10 5 pho- 
ton s _1 cm -2 sr _1 in the redshift interval 2 ^ z ^ 5 for 
2" pixels (see Fig. [6]). The strongest IGM emission comes 
from Hi Lya, provided we include emission from gas with 
10~ 3 cm~ 3 < n H < 10 _1 cm _1 . The highest surface bright- 
ness metal-line emission comes from C in and is about an 
order of magnitude fainter. The highest surface brightness 
regions for Civ, Sim, Siiv and Ovi lines are fainter than 
the brightest Cm by factors of a few. The Nv and Nevm 
lines, as well as Hen Ha, are substantially fainter but their 
maximum surface brightnesses still exceed 10 2 photon s _1 
cm -2 sr _1 for z — 2 (and 2" pixels). 

(ii) Lines from different ions trace gas in different tem- 
perature regimes. The gas responsible for the highest sur- 
face brightness metal-line emission typically has tempera- 
tures close to the peak of the corresponding emissivity curve. 
Hi Lya, Hen Lya, Cm, Sim and Siiv lines prevalently 
trace gas with T < 10 5 K near galaxies, making them excel- 
lent probes of cold accretion flows and the colder parts of 
galactic winds. Civ (T ~ 10 5 K), Nv (T ~ 10 51 K), Ovi 
(T ~ 10 5 ' 4 K), and Nevm (T ~ 10 5 ' 8 K) trace the warmer 
gas in haloes and filaments and thus constitute powerful 
tools to study the WHIM and the warmer parts of outflows. 

(iii) While the brightest Hi Lya emission arises exclu- 
sively in highly overdense gas (p/ (p) > 10 3 ), at least if we 
assume the gas to be optically thin, a much wider range of 
overdensities contributes to the brightest pixels of He n Lya 
and high-ionisation metal lines such as C iv, O vi and Ne vm. 
The spatial distribution of the emission is related to the ion- 
isation state of the emitting gas. Lower ionisation gas pro- 
duces clumpier and more clustered emission, while higher 
ionisation gas produces smoother and more diffuse emission 
profiles. The degree of ionisation decreases with increasing 
density but increases with the temperature. As the emis- 
sivity increases with the gas density, emission from lower- 
ionisation gas tends to be brighter. Indeed, we find for ex- 
ample that the beryllium-like atoms N iv and O v (which 
were not shown here because they only shift into the optical 
waveband for very high redshifts) are stronger than their 
higher-ionisation lithium-like counterparts. 

(iv) We have compared our fiducial simulation with mod- 
els that ignore SN feedback and/or metal-line cooling, with 
a model that includes highly mass-loaded "momentum- 



driven" galactic winds, and with a simulation that takes 
feedback from AGN into account. While the faint metal- line 
emission is very sensitive to the inclusion and implemen- 
tation of galactic winds (because these are responsible for 
enriching the low-density IGM with metals, where the weak 
emission is produced), the results for bright metal- line emis- 
sion, and even more so for both faint and bright Hi Lya 
emission, are remarkably robust to changes in these physi- 
cal processes. Accurate cooling rates are required to predict 
the bright metal-line emission with a precision better than 
a factor of a few, particularly for lines that peak at tem- 
peratures for which the cooling time is short, such as Ovi. 
AGN feedback primarily affects the emission from the cores 
of haloes, reducing the brightest metal-line emission by up 
to factors of a few. Highly mass-loaded galactic winds, on 
the other hand, can boost the brightest metal-line emission 
by similar factors. 

(v) According to our predictions, MUSE and KCWI 
should be able to detect H i Lya emission from regions with 
p > 10 3 pb up to at least z — 5 (provided our optically thin 
calculation is appropriate in this regime) , while H i Lya from 
less dense regions and the He n Ha line should be visible at 
lower redshifts. C iv, O vi and Si iv will all be detectable up 
to z ss 3, while Cm and Sim should be detectable up to 
z fa 4. The Nv and Nevm lines should also be within the 
capabilities of these instruments. We note, however, that 
lines with rest-frame wavelengths shorter than that of Lya 
(i.e. Cm, O vi, Ne vm, and Si in) may suffer substantial fore- 
ground absorption, which we have not accounted for here. 
O v lines at z as 2.2 and C in lines at z ~ 1 are roughly within 
the sensitivity of the FIREBALL experiment, and should be 
considered as possible contaminants of Lya emission. 

We conclude that several rest-frame UV emission lines 
from the high-redshift IGM will become detectable in the 
near future, possibly starting with the CWI which is already 
operating on Palomar. Observations of these lines have the 
potential to revolutionise studies of the interactions between 
high-redshift galaxies and their environment. Lower ioni- 
sation lines provide us with tools to image cold accretion 
flows as well as cold, outflowing clouds. Higher ionisation 
lines, on the other hand, will allow us to image the metal- 
enriched WHIM and galactic winds. Such observations, par- 
ticularly the simultaneous detection of multiple lines, will 
provide strong constraints on the density, temperature, and 
chemical composition of the gas. 
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APPENDIX A: ANGULAR RESOLUTION 

In this Appendix we test the effect of varying the angular 
resolution used to make the emission maps. Fig. I All com- 
pares the PDFs of the surface brightness of, from left to 
right, H i Lya, C iv, and O vi at z = 2 for maps with resolu- 
tions = 0.5" , 1" , 2" , 4" , 8" , 16" , and 32" , which correspond 
to physical scales varying from 3.1 kpc to 199 h~ x kpc. 

The distributions converge well for all but the high- 
est surface brightness values. The increase in the maximum 
brightness with decreasing pixel size reflects the fact that 
large pixel sizes smooth the strong emission coming from 
structures with an angular size smaller than that used to ob- 
serve them. The differences between the PDFs correspond- 
ing to different pixel sizes are much smaller for Ovi than 
for Civ and Hi. This is because Ovi emission arises in gas 
of lower densities, which will have larger coherence lengths 
(becaus e the Jeans s cale increases with decreasing density; 
see e.g. lSchavdkoOlT ). 



APPENDIX B: CONVERGENCE TESTS 

In this Appendix we present a number of tests aimed to 
verify how resolution affects the numerical predictions. We 
consider the effect of varying the numerical resolution in 
Appendix IB1I the box size in Appendix IB2I and the slice 
thickness in Appendix lB3l All tests use the REF model and 
results are shown only for the Ovi 1032 A line. 



Bl Numerical resolution 

To test the effect of varying the numerical resolution, we 
compare emission maps for three different realisations of the 
same simulation with particle numbers iV = 2x 128 3 , 2x256 3 
and 2x 512 3 . The corresponding dark matter particle masses 
are m DM = 4.1 x 10 s , 5.1 x 10 7 , and 6.3 x 10 6 h" 1 Mq, 
respectively. 

Fig. IB II shows the PDFs of the z = 2 surface brightness 
for, from left to right, H i Lya, C iv, and O vi for maps with 
2" angular resolution. All distributions are similar at high 
and intermediate surface brightnesses. Only the lowest res- 
olution run with N = 2 x 128 3 deviates strongly from the 
other two and only for < 10 -5 photon s _1 cm~ 2 sr _1 . We 
conclude that our conclusions are not significantly affected 
by the finite resolution of our simulations. 



B2 Box size 

We test the effect of varying the size of the simulated box, 
while keeping the mass resolution constant. The runs neces- 
sarily assume different initial conditions and thus allow us to 
test how the large-scale structure affects the flux statistics. 

In Fig. IB2I we compare results for H i Lya, C iv, and 
Ovi at z = 2 using simulations with dark matter particle 
mass rriDM = 6.3 x 10 6 /i -1 Mq and box sizes of 25, 12.5, and 
6.25 ft -1 Mpc (comoving), respectively. The maps used to 
calculate the PDF have 2" angular resolution and 6.25 h~ l 
Mpc thickness, consistent with the box size of the smallest 
simulation. The results are converged with respect to the 
size of the simulation volume, although larger boxes sample 
the flux PDF to higher values. 



B3 Slice thickness 

In this Section we describe how varying the thickness of the 
slices we cut through the simulation box affects the shape 
of the flux PDF. In Fig. IB3I we compare the z — 2 surface 
brightness PDFs of, from left to right, Hi Lya, Civ, and 
Ovi for three different, evenly-spaced slice thicknesses: 25, 
5 and 1 h~ 1 Mpc. To facilitate the comparison, we have 
normalised the surface brightness PDF by the slice thick- 
ness Ax. Throughout this work, we have assumed a slice 
thickness of 25 h' 1 Mpc. 

Except for C iv, the results appear to be independent 
of the slice thickness (within limits and in the detectable 
regime) when the PDF is normalised by the slice thickness. 
In other words, the normalisation of the PDF is propor- 
tional to the slice thickness. This happens when the slice 
is thick enough to fully contain collapsed objects, but suf- 
ficiently small to avoid the superposition of multiple struc- 
tures, which would strongly increase the flux in those pixels 
where more than one structure could be found along the line 
of sight. For Civ there is a small difference between 5 and 
25 ft -1 Mpc, but not between 1 and 5 h^ 1 Mpc, indicating 
that projection effects slightly enhance the flux for 25 h~ l 
Mpc slices. 

The shape of the PDF at the faint end, on the other 
hand, depends on the level of the flux in low-density regions 
and on the number of low-density structures in each slice, 
which is dependent on the slice thickness. As a consequence, 
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Figure Al. As Fig. [5] but for Hi Lya (left panel), Civ (middle panel), and Ovi (right panel) emission as a function of the angular 
resolution of the emission maps. The pixel sizes are i9 = 0.5", 1", 2", 4", 8", 16", and 32", corresponding to physical sizes of 3.1, 6.2, 
12.4, 24.8, 49.7, 99.3, and 199 h' 1 kpc, respectively, at z = 2. The maximum predicted surface brightness increases steadily with the 
angular resolution. The results are less sensitive to the pixel size for lines that arise in lower density gas, such as O VI. 
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Figure Bl. As Fig. [5] but for Hi Lya (left panel), Civ (middle panel), and Ovi (right panel) emission as a function of the numerical 
resolution of the simulation. The particle mass increases by a factor of eight and the force resolution decreases by a factor of two when 
the number of particles decreases by a factor of 2 3 . Our fiducial run is converged with respect to the numerical resolution. 
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Figure B2. As Fig. [5] but for Hi Lya (left panel), Civ (middle panel), and Ovi (right panel) emission for simulations with box sizes 
of 25, 12.5, and 6.25 h~ 1 Mpc (comoving) and a fixed numerical resolution. Different from Fig. [5] all maps are for 6.25 h~ 1 Mpc thick 
slices. The results are converged with respect to the size of the simulation volume. 
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Figure B3. The surface brightness PDFs for Hi Lyo (left panel), Civ (middle panel), and Ovi (right panel), normalised by the slice 
thickness Ax, for maps with varying slice thickness: 25, 5 and 1 /i" 1 Mpc. The number of slices used to calculate each PDF is 1, 5 and 
25 respectively. The test is for the REF model with box size 25 h~ l Mpc at z = 2. All maps assume the same angular resolution of 2". 
For Hi Lyo and O VI the PDF is proportional to the slice thickness in the high flux regime. This indicates that the objects that produce 
the highest fluxes in the distribution are fully contained within a slice. For Civ there is a small difference between 5 and 25 /i -1 Mpc, 
but not between 1 and 5 h~ 1 Mpc, indicating that projection effects slightly enhance the flux for 25 h^ 1 Mpc slices. 

the fraction of pixels at the lowest fluxes increases for de- 
creasing slice thickness. 

This paper has been typeset from a Tjj)X/ I^TjTjX file prepared 
by the author. 
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